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Abstract 

In this paper we investigate the quantum and classical dynamics of a single 
trapped ion subject to non-linear kicks derived from a periodic sequence of 
Gaussian laser pulses. We show that the classical system exhibits diffusive 
growth in the energy, or 'heating', while quantum mechanics suppresses this 
heating. This system may be realised in current single trapped-ion experi- 
ments with the addition of near-field optics to introduce tightly focussed laser 
pulses into the trap. 
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I. INTRODUCTION 



Recent experiments on the nonlinear dynamics of cold trapped atoms have provided a 
remarkable verification of key theoretical ideas in the subject of quantum chaos includ- 
ing dynamical localisation and the effect of decoherence in restabilising the classical limit. 
In all these experiments however the observed results are obtained from a large ensemble 
of single atom experiments which run in parallel but independently. Until now, there have 
been no experiments which investigate quantum chaos of a single particle monitored over a 
period of time. In contrast, in the related field of ion-trapping, technological advances now 
enable a single ion to be trapped, cooled to the ground state of the trap, and monitored, 
almost without error || . The quantum dynamics of the centre of mass motion of the ion is 
extremely well described by a three dimensional harmonic potential. In some experiments, 
two degrees of freedom are very tightly bound and the interesting harmonic motion takes 
place in a single degree of freedom. Of course this system is integrable. However if this 
degree of freedom is subject to a periodic nonlinear potential, chaos may result. In this 
paper we investigate the quantum and classical dynamics of a single trapped ion subject to 
non-linear kicks derived from a periodic sequence of Gaussian laser pulses. This system may 
be realised in current single trapped-ion experiments with the addition of near-field optics 
to introduce tightly focussed laser pulses into the trap. Another suggestion for investigating 
quantum chaos in a single trapped ion has recently been suggested by Berman et al J4j. 

The recently achieved ability to engineer dynamics for a single trapped ion has followed 
from the potential application of this system for quantum computational gates. As such 
these systems necessarily operate at the quantum level and provide the ideal experimental 
context to test quantum nonlinear dynamics. Indeed such experiments will ultimately involve 
the ability to follow the quantum dynamics of many trapped ions with complex many- 
body interactions introduced by externally imposed time dependant Hamiltonians. We 
thus believe it timely to consider tests of quantum chaos which can be made with current 
technology. 

In section II we define the classical dynamical system and give a detailed analysis of the 
classical motion and the transition to chaos. In section III we give a quantum description of 
the problem and in section IV show that this system exhibits a suppression in the diffusion 
of momentum and position. In other words the total energy of the trapped ion is localised 
in contrast to classical diffusive heating of the motion. Finally in Section V we discuss a 
possible physical realisation of the system. 



II. CLASSICAL MAP 



The system has the Hamiltonian 

H = f - + ^-x 2 + Ke^T T 8(t - nT), 

2171 2 n^oo 

where m is the mass of the ion trapped in a harmonic potential of frequency Q. The ion 

~ 2 

is subject to a periodic sequence of laser pulses with period T. The ne~ ax term in the 
Hamiltonian describes the potential felt by the ion due to the Gaussian structure of the 
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laser. If we rescale time, position and momentum by letting i = tT, x = -^x and p = ^=jo 
the Hamiltonian rescales to 

rji OO 

H(x,p,t) = ^-H(5;,p,t) = ^(p 2 + x 2 )+ke-* 2 £ 5(t-n). (1) 

Tl— — OO 

where co = QT and k = are dimensionless parameters. The new variables x, p and t 
together with the Hamiltonian, H, are also dimensionless. 

Between kicks the system has the solution (x,p) = (xq cos cot + po sin out, — Xq sin cot + 
Po cos cot). The effect of the kick is to add a position dependent shift in the momentum of 
2kxe~ x . Denoting the mapping by F we can write it as the composition of a kick, K, and 
a linear rotation, W, 

F = Wo K 

where, 

2 

K(x,p) = (x, p + 2kxe~ x ) 

W(x, p) = (x cos co + p sin co, —x sin co + p cos co) 

Hence F maps (x,p) from just before a kick to one period later. 
The fixed points of F are at 



(x,p) = (0,0), (^±^ln(A;cot|),Ttan|yin(A;cot|)^ . 

The origin is stable if k cot ^ < 1 and —A; tan f < 1. When k cot ^ > 1 it becomes unstable 
via a pitchfork bifurcation which creates the second two fixed points. These are stable 
for A; cot ^ < exp(|csc 2 ^). When — fctan^ > 1 the origin becomes unstable via a period 
doubling bifurcation and two period 2 orbits are created at 



(x, p) = (±^ln(-fctan^), =F cot ^ ln(-k tan^ . 

These are stable for —A; tan ^ < exp(^ sec 2 ^). 

One can also find two sets of period 4 orbits which exist for —A; tana; > 1. The first set 
is at 

(x, p) = (±x, ±X cot |) , (±x, TX tan |) (2) 



where x — y hi(— A; tan a;). They are stable if —ktanu < exp(|| seco;|). A second set lie in 
between these orbits at 

(x,p) = (±x,±xcotco),(0,±xcscco). 

These are always unstable. 

It is instructive to study the process of creation and destruction of these period 4 orbits 
in more detail. In Fig. 1 we have drawn a bifurcation diagram. The period 4 orbits exist 
for parameter values lying in the shaded regions. The regions of lighter shade show where 
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the first set of orbits (0) are unstable. In figures 2-6 a sequence of phase space pictures are 
drawn for values of uj along the line k = 2. Fig. 2 is for ^- = 0.24 and shows typical phase 
space structure for the system. As uj is increased, the arms of the star-shaped chaotic region 
grow in size. This can be seen in Fig. 3 where ^- = 0.248. At ^- = \ (Fig. 4) these arms 
are infinitely long chaotic channels which divide the phase space into four regions, creating 
the period 4 orbits at infinity. Fig. 5 shows the period 4 orbits just after they are created 
(^- = 0.252). They now move towards the origin as uj is increased further. On this journey 
the first set (|2]) shed their stability via a pitchfork bifurcation but then regain it before 
destroying at the origin. Fig. 6 shows the orbits just before destruction (j- = 0.42). 

It is simple to see how the period 4 orbits were created at infinity when uj — The kick 
has little influence on orbits here and thus F reduces to simple linear rotation with period 
4. In general, it can be shown that an orbit of period p is created at infinity when 

2nq 
UJ = , 

P 

and then destroyed at the origin when 

cos uj + k sin uj = cos , 

P 

where q and p are natural numbers with a greatest common divisor of one. The condition 
for destruction is found by looking at the eigenvalues of the tangent map at the origin and 
then equating this linear rotation with the period of the orbit. 



III. QUANTUM MAP 

To construct the quantum map we start with the rescaled Hamiltonian (P and define a 
dimensionless Planck's constant, k, via the commutation relation for position and momen- 
tum 

x.p\ = [x,p\ = i = IK. 

muj muj 

The time evolution of an initial state, \ip n ), from just before a kick through to one period 
later is given by 

= ex? J\{t)dt)\r) (3) 

= exp (p 2 + x 2 )) exp \V) (4) 

= Pm (5) 

Hence the Floquet operator, F, defines the quantum map. Now, defining the annihilation 
and creation operators to be 

x + ip) and a) = (x — ip) 



2k V2k 
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respectively, and using the simple harmonic oscillator eigenstates 

\n) = -^=a ]n \$) n = 0,l,2... 

vn! 



as an orthogonal basis we can rewrite equation as 

n+ i „ 

Sn ± mk^k 

where cj^ = (m\ifj n ) and 

i^nm = (n|F|m) 



-icj(n+i) 



— i/c 



,.2 



3^(n| exp £~ e 1 J l m )- 

The last term is found by taking the exponential of the matrix with components 
/ ,—ik _ x 2, v —ik f°° „ / x \ TT ( x \ -(i+^)^ 2 



n 



k ' k 3/2 y/^n\m\2 n+m J-oo n \Vk J m \Vk 



m) = „,„ - ; ; / H n [-= )H m [-= e rfx 



— A;(n + ?n — l)(n + m — 3) ■ ■ • 1 / — \ 2 / 1 — n — ml 1 



for n + m even and vanishing otherwise. Here if n are Hermite polynomials and 2F1 is the 
hypergeometric function. Note that F nm = whenever n + m is odd. This means that even 
and odd parity states do not couple under F and thus evolve independently. 



IV. LOCALISATION 

We now show numerically the presence of dynamical localisation in the system. Or, 
more precisely, we show that classical diffusion is suppressed when the system is evolved 
quantum mechanically. For this, we have chosen uj = ir(3 — y/E) and k = 8. In this 
parameter regime a large chaotic sea centred at the origin consumes the phase space (see 
Fig. 7). The initial state was chosen to be |0), which has a Husimi probability density of 

|<*|o>| 3 = e -w a (6) 

in phase space. Here z = —}=(x + ip) and \z) are the coherent states defined as 

\z) = e~2N ^ —i=\ri). 

Thus our initial state is a highly localised Gaussian hump centered at the origin. This 
was then evolved forward using 1800 of the even basis states. Fig. 8 shows the average 
dimensionless energy, (x 2 +p 2 ), after each kick. The light gray is for k = 0.5 and the dark 
gray is for k = 0.2. The energy under classical evolution is shown in black. Here an initial 
density (||) with k = 0.2 was chosen. One can clearly see in this figure that diffusion is 
suppressed after about 100 kicks. The procedure was repeated using only 1000 even basis 
states to confirm accuracy and it was found that the difference in the energies did not exceed 
10 -8 until after 1000 kicks. Thus the localisation is truly a property of quantum mechanics. 
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V. DISCUSSION AND CONCLUSION 



What are the physical requirements to realise this system in current ion trap experiments? 
Consider a 9 Be + ion such as used in the NIST experiments ||, with a harmonic frequency 
in the relevant direction of Q = 1 MHz. The key parameter which determines the effective 
Planck constant is the parameter a. If one uses a focused laser beam, a typical value for 
this parameter is 10 10 m 2 and a resulting effective Planck's constant of k = 6.7 x 1CT 5 , 
which is hopelessly too small. On the other hand if we use a near-field probe, as used in 
near-field optical scanning microscopy (NOSM), to inject the field we can get a value as high 
as a = 10 14 m 2 with a typical probe tip of diameter 10 nm. This corresponds to an effective 
Planck's constant for 9 Be + of about k = 0.7, which is more promising. To achieve a value 
for the kick parameter of the order used above we would need to focus a few nanowatts into 
the NOSM probe which is quite typical. This would correspond to an intensity of about 
1000 mW cm -2 at the ion. If we choose the kick period to be of the order of 10 ms, the 
kick parameter, k, has a value of the order of unity. We conclude that this experiment is 
possible for a current single trapped-ion experiment with the addition of near-field optical 
fibre probes. 

The next question we need to ask of such a system is how are we to observe the motion 
of the ion? Fortunately the current single trapped-ion experiments are designed precisely 
to enable careful monitoring of the motion states. The details are described in reference [|J. 
The basic idea is to map the motion states onto particular internal states of the ion which 
are then probed by a fluorescent shelving technique. In particular it is possible to measure 
the centre-of-mass energy of the ion in the trap. Each measurement however destroys the 
quantum state of the ion at that time, so repreparation of the ion initial state is required. 
One then needs to perform repeated experiments for differing number of kicks before reading 
out the centre-of-mass energy. In this way it is possible to monitor the energy of the motion 
as a function of kick number. Dynamical localisation of the motion energy of the ion could 
thus be observed. 

Finally we need to ask if it is feasible to prepare the initial states we have used in this 
paper. Again reference || shows that it is possible to prepare the ion in the ground state 
of the harmonic trap, so this part is relatively easy. Laser pulses may then be used to 
displace this minimum uncertainty state anywhere in the phase plane. This ability to place 
a localised state anywhere in the phase plane would enable a detailed study of mixed chaotic 
and regular phase space structures. Unfortunately in the current experiment unwanted stray 
linear potentials cause a heating of the ion and thus it does not stay in the ground state for 
long, but rather undergoes a diffusive motion in the phase plane |||7|]. A very considerable 
amount of effort is currently being devoted to removing this unwanted heating so that 
trapped ions can be used in a quantum logic gate. We thus expect this problem to be solved 
or at least significantly mitigated. 

Needless to say this is not an easy experiment. Introducing the near field probe close 
to the ion will cause additional unwanted van der Walls forces to be exerted on the ion. 
However these forces, while making a detailed comparison to experiment more difficult, will 
not effect the generic transition to chaos described above so long as they remain weak. The 
heating of the ion due to stray linear potentials will remain a problem to some extent. Such 
fluctuating forces are a source of decoherence and thus will tend to destroy localisation. 
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Taking a longer view however the ease with which decoherence can be induced via this 
mechanism should enable a detailed study of the effect of noise on dynamical localisation to 
be made, thus turning a bug into a feature. 



VI. FIGURE CAPTIONS 



FIG. 1. Bifurcation diagram for the period 4 orbits (0). The stable region is in dark 
and the unstable in light. 

FIG. 2. Phase portrait for k = 2 and £ = 0.24. 

FIG. 3. As for Fig. 2 except ^ = 0.248. 

FIG. 4. As for Fig. 2 except £ = \. 

FIG. 5. As for Fig. 2 except £ = 0.252. 

FIG. 6. As for Fig. 2 except £ = 0.42. 

FIG. 7. Phase portrait for k = 8 and to = 7r(3 — Vo). 

FIG. 8. The average dimensionless energy after each kick for the initial state |0) with 
k = 0.5 (light gray), k = 0.2 (dark gray) and classical evolution (black). 
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